Astronomy & Astrophysics manuscript no. aa20377 


©ESO 2012 


October 15, 2012 





Letter to the Editor 



Four temporary Neptune co-orbitals: (148975) 2001 XA 25 5, (310071) 
2010 KR 59 , (316179) 2010 EN 65 and 2012 GX 17 * 

C. de la Fuente Marcos and R. de la Fuente Marcos 



(N 



Universidad Complutense de Madrid, Ciudad Universitaria, E-28040 Madrid, Spain 
Received 13 September 2012 / Accepted 10 October 2012 

ABSTRACT 



o 



Context. Numerical simulations suggest that Neptune primordial co-orbitals may significantly outnumber the equivalent population 
1 hosted by Jupiter, yet the objects remain elusive. Since the first discovery in 2001 just ten minor planets, including nine Trojans and 

one quasi-satellite, have been positively identified as Neptune co-orbitals. In contrast, Minor Planet Center (MPC) data indicate that 
more than 5,000 objects are confirmed Jupiter co-orbitals. On the other hand, some simulations predict that a negligible fraction of 
passing bodies are captured into the 1 : 1 commensurability with Neptune today. 

Aims. Hundreds of objects have been discovered in the outer solar system during the various wide-field surveys carried out during 
| the past decade, and many of them have been classified using cuts in the pericentre and other orbital elements. This leads to possible 

misclassifications of resonant objects. Here, we explore this possibility to uncover neglected Neptune co-orbitals. 
Methods. Using numerical analysis techniques, we singled out eleven candidates and used A'-body calculations to either confirm or 
reject their co-orbital nature. 

Ql^ ' Results. We confirm that four objects previously classified as Centaurs by the MPC currently are temporary Neptune co-orbitals. 

, (148975) 2001 XA255 is the most dynamically unstable of the four. It appears to be a relatively recent (50 kyr) visitor from the 

scattered disk on its way to the inner solar system. (310071) 2010 KR59 is following a complex horseshoe orbit, (316179) 2010 EN65 
i-C 1 is in the process of switching from L 4 to L 5 Trojan, and 2012 GX I7 is a promising L 5 Trojan candidate in urgent need of follow-up. 

The four objects move in highly inclined orbits and have high eccentricities. These dynamically hot objects are not primordial 1:1 
librators, but are captured and likely originated from beyond Neptune, having entered the region of the giant planets relatively recently. 
Conclusions. Casting doubt over claims by other authors, our results show that Neptune can still efficiently capture co-orbitals for 
short periods of time and that the cuts in the orbital elements are unreliable criteria to classify objects orbiting in the outer solar 
system. As in the case of Jupiter Trojans, our results suggests that Neptune's L5 point is less stable than L4, in this case perhaps due 
to the influence of Pluto. 
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C*") [ 1. Introduction Confirmed Jupiter co-orbitals currently amount to more than 

(— ^ . 5,000 objectsj; theoretical and numerical expectations predict 

t-j . Numerical simulations predict that Neptune may have retained eyen , numbers for Ne ^ ^ fa ± ^ harder tQ 

■ a significant amount of its primordial co-orbital minor planet detect becaufje ^ m farther Qn ^ Qther hand; ft ^ 

t-H ■ P°P" latl0 f' T«r ng Ja f ™ q ^ a S1 " Sat o S J> ol? commonly thought that Neptune cannot efficiently capture ob- 
. Wisdom 1993; Wiegert et al 2000; Nesvorny & Dones 2002; & commen L ability even for J n p l liods of 

Marzan et al. 2003; Kortenkamp et al. 2004; Dvorak et al. time (Horner & Evans 2006) 

k> 2007; Nesvorny & Vokrouhlicky 2009; Zhou et al. 2009, 2011; ' . t , . , , . , B , , ftU t , 
X T , , . , onnri ~ nin onl ,x r™. c . * T . L ■ During the past decade, wide-field surveys of the outer solar 
•\> Lykawka et al. 2009, 2010, 2011). The first Neptune Troian, ^ , ■ ■, 
H, ... , , , „, . J system have found hundreds of objects passing in the neighbour- 
ed 2001 QR322, was serendipitously discovered by Chiang et , , j. XT r , , .„ , . 1 • ^ 
1 al. (2003). 2004 UP 10 , 2005 TN53 and 2005 TO 74 followed hood of Neptune, many of which were classified using cuts in the 

(Sheppard & Trujillo 2006). Then 2006 RJ103 and 2007 VL 305 , 
all of them leading 60° ahead of Neptune (i.e. they are L 4 



pericentre, q, and other orbital elements (for example, the MPC 
defines that Centaurs must have a perihelion larger than Jupiter's 



„ . . „, . „ 6 , „ T „ . F ,„„„ , _ j- j orbit and a semi-maior axis shorter than Neptune s). This leads 

Troians). Shortly after, the first L5 Troian, 2008 LCis, was found , .„ c t u- « /ou 1 \ n ,,v y . 

to misclassifications of resonant objects (Shankman 2012). It is 



(Sheppard & Trujillo 2010), followed by the second one, 2004 

KV\«. These eight Trojans have recently been joined by one ,,,,,, ,- , •,- , 

, ° . „ ...i /-, nm im inna mu /j i t- identified and that they are, in tact, trapped (even it temporarily) 

short-term quasi-satellite, (309239) 2007 RWin (de la Fuente . , V1 _ XT „ T i 1 



possible that some of these objects may have not been properly 



Marcos & de la Fuente Marcos 2012a) and yet another L5 
Trojan, 2011 HMi 02 , (Tholen et al. 2012; Parker et al. 2012). 



in a 1 : 1 commensurability with Neptune. Here, we explore this 
possibility and try to uncover misidentified Neptune co-orbitals. 
In this letter, we use A^-body simulations to confirm the possible 
Send offprint requests to: C. de la Fuente Marcos, e-mail: co-orbital nature with Neptune of a number of objects currently 
nbplanet@fis.ucm.es classified as Centaurs by the Minor Planet Center (MPC). The 



Tables l2l4l and Figures l4l6l are only available in electronic form via 



http://www.edpsciences.org http://www.minorplanetcenter.net/iau/lists/JupiterTrojans.html 
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numerical model is described in the next section and the results 
are presented in Sections 3 to 6. The results are discussed and 
our conclusions summarized in Section 7. 

2. Candidates and numerical experiments 

Following Mikkola et al. (2006) and to study the librational 
properties of any possible candidates, we define the relative de- 
viation of the semi-major axis from that of Neptune by a — 
(a-aN)/apf, where a and are the semi-major axes of the object 
and Neptune, respectively, and also the relative mean longitude 
A r = A - An, where A and A^ are the mean longitudes of the ob- 
ject and Neptune, respectively. If A r oscillates around 0°, we call 
the object a quasi-satellite; a Trojan is characterized by A r os- 
cillating around +60° (L 4 Trojan) or -60° (or 300°, L 5 Trojan); 
finally, an object oscillating with amplitude > 180° follows a 
horseshoe orbit in a frame of reference rotating with Neptune 
(see, e.g., Murray & Dermott 1999). Objects that switch libra- 
tion from the Lagrangian point L4 to L5 (or vice versa) are called 
jumping Trojans (Tsiganis et al. 2000). Searching for Neptune 
co-orbital candidates among known minor bodies with relative 
semi-major axis \a - a^\ < 2 AU, we found 21 objects, 10 of 
which were already documented as co-orbitals. We regarded the 
remaining 1 1 as candidates and then used /V-body calculations 
to either confirm or reject their putative co-orbital nature. The 
orbit computations were completed using a Hermite integration 
scheme (Makino 1991; Aarseth 2003). The /V-body code is pub- 
licly available from the IoA web sitq3- Our calculations include 
the perturbations by the eight major planets, treating the Earth 
and the Moon as two separate point masses, the barycentre of 
the Pluto-Charon system, and the three largest asteroids. For ac- 
curate initial positions and velocities we used the heliocentric 
ecliptic Keplerian elements provided by the JPL on-line solar 
system data servic^l (Giorgini et al. 1996) and initial positions 
and velocities based on the DE405 planetary orbital ephemerides 
(Standish 1998) referred to the barycentre of the solar system. In 
addition to the orbital calculations completed using the nominal 
elements in Table Q] we have performed 50 control simulations 
(for each object) with sets of orbital elements obtained from the 
nominal ones and the quoted uncertainties (3cr) when available. 
Additional details can be found in de la Fuente Marcos & de la 
Fuente Marcos (2012b). We have also validated our simulations 
against previous work by integrating the orbits of the known nine 
Neptune Trojans (orbital elements in Tables |2]|4|, see Figs. |4jj6j 
results are consistent with those from other authors. 
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Fig. 1. Motion of the objects in a coordinate system rotating with 
Neptune. The orbit and position of Neptune are also plotted. The 
time intervals displayed are 2001 XA 25 s (-1.5, 1.5) kyr, 2010 
KR 5 y (5, 20) kyr, 2010 EN 65 (-20, 20) kyr, and 2012 GX i7 (-50, 
50) kyr. 
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Fig. 2. Resonant evolution. The a parameter (see text for details) 
as a function of the relative mean longitude during the time in- 
terval (-50, 50) kyr, except for 2012 GX 17 , (-100, 100) kyr. 
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3. (148975) 2001 XA 255 : very dynamically unstable 

(148975) 2001 XA 255 was discovered on December 9, 2001 at 
Mauna Kea Observatory by Jewitt et al. (2002). It moves in a 
very eccentric orbit (0.68) with perihelion just inside the orbit 
of Saturn and significant inclination (12.6°). Its aphelion is in 
the trans-Neptunian belt. It is classified as an inactive Centaur 
(Jewitt 2009). Its colours are neutral (Fraser & Brown 2012). 
(148975) 2001 XA255 is the most dynamically unstable of the 
four objects studied here with an e-folding time (or character- 
istic timescale on which two arbitrarily close orbits diverge ex- 
ponentially) of about 300 yr. It is perturbed by Saturn, Uranus, 
and Neptune. It appears to be a relatively recent (50 kyr) visi- 
tor from the scattered disk on its way to the inner solar system. 
Calculations suggest that it came from beyond 100-200 AU and 

2 http://www.ast.cam.ac.uk/~sverre/web/pages/nbody.htm 

3 http://ssd.jpl.nasa.gov/7planet_pos 



it is currently a very short-lived horseshoe librator (see Figs. Q] 
and|2]i. It entered the horseshoe orbital path about 10 kyr ago and 
will leave its present dynamical state 2 kyr from now (see Fig. 
0; very close encounters with Neptune are possible. It has pre- 
viously been identified by Gallardo (2006) as a transition object 
affected by a 1:1 mean motion resonance with Neptune and 1:2 
with Uranus and was classified by Bailey & Malhotra (2009) as 
a difussing Centaur. 

4. (310071) 2010 KR 59 : a transient horseshoe 

Discovered on May 18, 2010 by NASA's Wide-field Infrared 
Survey Explorer (WISE) (Scotti et al. 2010). The WISE tele- 
scope scanned the entire sky in infrared light from January 
2010 to February 2011 (Wright et al. 2010). (310071) 2010 
KR59 moves in a very eccentric orbit (0.57) with perihelion out- 
side Saturn's orbit and significant inclination (19.67°). Little is 
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Table 1. Heliocentric Keplerian orbital elements of the objects studied in this research. Values include the 1-cr uncertainty when 
available. (Epoch = JD2456200.5, 2012-Sep-30.0 except for 2012 GX n that is JD2456020.5, 2012-Apr-3.0; J2000.0 ecliptic and 
equinox. Sources: JPL Small-Body Database and AstDyS-2.) 





(148975) 2001 XA 255 


(310071) 2010 KR59 


(316179) 2010 EN 65 


2012 GX 17 


semi-major axis, a (AU) = 
eccentricity, e = 
inclination, < (°) = 
longitude of the ascending node, fi (°) = 
argument of perihelion, at (°) = 
mean anomaly, M (°) = 
perihelion distance, q (AU) = 
absolute magnitude, H (mag) = 


28.8587±0.001() 
0.676579±0.000010 
12.62072±0.00002 
105.92560±0.00014 
90.3433±0.0004 
5.3117±0.0002 
9.33350±0.00010 
11.2 


29.970±0.004 
0.56676±0.00006 
19.67387±0.00005 

46.7483±0.0003 
108.2697±0.0010 

5.6713±0.0014 

12.9842±0.0008 
7.7 


30.719±0.002 
0.31492±0.00004 
19.24716±0.00003 
234.2657±0.0003 
225.236±0.003 
36.831±0.005 
21.0450±0.0009 
6.9 


30.1330172 
0.4149134 

35.29849 
207.75190 

94.00361 
298.4401 1 
17.6304246 
7.8 
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Fig. 3. Time evolution of various parameters for 2001 XA255, 2010 KRg9, 2010 EN65, and 2012 GX17. The distance to the object 
from Neptune, diy, is shown in panel A and the value of the Hill sphere radius of Neptune, 0.769 AU, is displayed. The resonant 
angle, A r , is shown in panel B for the nominal orbit in Table [TJ (thick red line) along with one of the control orbits (thin blue line). 
This particular control orbit was chosen to lie close to the 3-<x limit so that its orbital elements are distinct from the nominal ones 
with the exception of 2012 GX17. For that object the control orbit elements deviate by 1 % from the nominal values in TableQ] which 
is at least 25 times higher than those of the other objects. The orbital elements are shown in panel C, a, with the current value of 
Neptune's semi-major axis, panel D, e, and panel E, i. 



known about this object with the exception of its absolute magni- 
tude of 7.7 (~22 apparent). This suggests a relatively large object 
with an estimated diameter of about 100 km. Its e-folding time is 
close to 10 kyr, and it is now pursuing a complicated horseshoe 
orbit (Fig. [2j in its way to become a short-term quasi-satellite as 
depicted in Fig.Q] It will leave the 1:1 resonance with Neptune 
in a few 100 kyr. Its current horseshoe state started about 100 kyr 
ago. Dangerously close encounters with Neptune are possible. 



5. (316179) 2010 EN 65 : a jumping Trojan 

(316179) 2010 EN 65 was discovered on March 7, 2010 by 
Rabinowitz and Tourtellotte using the 1 .3-m reflector from Cerro 
Tololo (Lowe et al. 2010). Upon discovery, the orbit of the ob- 
ject was described as Neptune-like, but its relative mean longi- 
tude with respect to Neptune was almost 180°. Up to 19 pre- 



covery images were found shortly after, the first one acquired 
on November 4, 1989 by the Digitized Sky Survey (DSS) from 
Palomar Mountain. It is the largest of the four objects studied 
here, perhaps with a diameter of nearly 200 km. It moves in an 
eccentric orbit (0.31) with perihelion just outside Uranus's orbit 
and significant inclination (19°). It was classified by Rabinowitz 
et al. (2012) as a scattered-disk object. Our calculations show 
that its e-folding time is close to 3 kyr and it is currently moving 
from the Lagrangian point L4 to L5 (see Figs. [TJ [2] and [3}. This 
result is very robust as all the control orbits follow the same or- 
bital behaviour within 40 kyr of the current epoch. We classify 
this object as a jumping Trojan. 
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6. 2012 GX| 7 : a promising U Trojan candidate 

The minorplanet 2012 GXn was discovered by Pan-STARRS 1, 
Haleakala, Hawaii on March 14, 2012 and was reobserved from 
Magdalena Ridge Observatory, Socorro, New Mexico on March 
18. As a recently discovered object and in sharp contrast with the 
previous three, its orbit is poorly known and it is mainly included 
here to encourage follow-up observations: its orbital solution is 
based on ten observations with an arc length of just four days. 
The physical properties of this object are unknown with the ex- 
ception of its absolute magnitude of 7.8 (~22 apparent). This 
suggests a relatively large object with an estimated diameter in 
the range 55-180 km (for an albedo range of 0.5-0.05 and an as- 
sumed density of 2000 kg irT 3 ). Its period, 165 yr, matches that 
of Neptune, 164.79 yr, therefore, it appears to follow a 1:1 reso- 
nant orbit with Neptune (a = 30. 13 AU), yet it has been classified 
as a Centaur by the MPC. Our calculations for the nominal orbit 
in Table Q] indicate that the value of the relative mean longitude 
of 2012 GX n librates around -60° with an amplitude of 40-55° 
and a period of about 10 kyr; therefore, it appears to be a trail- 
ing Trojan (see Figs.[T]|2]and[3]). Although its e-folding time is 
nearly 10 kyr, it has been an L5 Trojan for tens of kyr and simu- 
lations suggest that it will be leaving its current dynamical state 
to become a horseshoe librator in about 50 kyr. 201 1 HM102 was 
the highest inclination Neptune Trojan known but this candidate 
has even higher orbital inclination (z = 35°). Due to its poorly 
known orbit, we must insist that the object is a mere Trojan can- 
didate (and in dire need of follow-up observations) although all 
studied control orbits (with errors below 1%) give consistent re- 
sults. Our calculations show that 2012 GXn may have been a 
passing minor body (perhaps from the scattered-disk) that even- 
tually evolved dynamically into a trailing Neptune Trojan. We 
also find that compared to Neptune's other Trojans, 2012 GXn 
is significantly more perturbed by Uranus. 

7. Discussion and conclusions 

Compared to the eight "classical" Neptune Trojans (2001 QR322, 
2004 UP 10 , 2005 TN 53 , 2005 TO 74 , 2006 RS m , 2007 VL305, 
2008 LCig, and 2004 KVi 8 ) and the new one (201 1 HM 102 ), the 
objects described here (even if they are currently co-orbitals of 
Neptune) are much more dynamically hot and their characteris- 
tic lifetime as co-orbitals is relatively short, but comparable to 
that of the most unstable of the other nine objects (2004 KVis, 
see Figs. |U|6]l. Dynamically speaking, they are similar to the re- 
cently identified temporary Neptunian quasi-satellite (309239) 
2007 RW 10 (de la Fuente Marcos & de la Fuente Marcos 2012a). 
The clue to the true origin of all these five objects (and, perhaps, 
2004 KVi 8 ) could be (148975) 2001 XA 255 , which clearly came 
from beyond 100-200 AU. This conclusion is similar to that in 
Horner & Lykawka (2012) for the case of the Neptune L5 Trojan 
2004 KVig. On the other hand, (316179) 2010EN 65 could be the 
brightest and 2012 GXn the highest inclination Neptune Trojan. 

Out of more than 5,000 confirmed Jupiter Trojans there 
is a well-documented asymmetry between the number of bod- 
ies in the leading and trailing populations: the Sloan Digital 
Sky Survey (SDSS) puts the L 4 :L 5 Trojan ratio at 1.6+0.1 
(Szabo et al. 2007), the Subaru telescope survey gives nearly 
1.8 (Nakamura & Yoshida 2008; Yoshida & Nakamura 2008) 
and the latest results from the WISE/NEOWISE mission give 
1.4+0.2 (Grav et al. 201 1). Prior to this research the L^.L^ ratio 
for Neptune was 6/3=2. Our calculations suggest that Neptune's 
L5 point is less stable than L4, apparently due to the influence of 
Pluto. All trailing Trojans undergo close approaches within 0.3 



AU of Pluto. Close encounters with the dwarf planet can send 
trailing Trojans into complex paths away from the usual tadpole 
orbit, inducing jumping Trojan events, temporary quasi-satellite 
episodes, or throwing the object on a horseshoe orbit. An asym- 
metry between the Neptunian L4 and the L5 swarms was found 
in numerical integrations by Holman & Wisdom (1993) but, at 
that time, it was not understood why that asymmetry existed. 

It is sometimes claimed that Neptune cannot currently effi- 
ciently trap objects in the 1:1 commensurability even for short 
periods of time. Our results argue that, contrary to this view, 
Neptune is still actively capturing temporary co-orbitals. This is 
consistent with recent findings by Kortenkamp & Joseph (201 1). 
These authors conclude that the nearly 2:1 mean motion reso- 
nance between Uranus and Neptune facilitates capture of new 
Neptune co-orbitals. Most of the objects described here are rel- 
atively large and a much more numerous population of smaller 
objects similar to them may exist. We also confirm that dynam- 
ical classifications based on cuts in the pericentre and other ele- 
ments are not very reliable in the case of resonant objects. 
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Table 2. Heliocentric Keplerian orbital elements of Neptune's L4 Trojans (I): 2001 QR322, 2004 UP10, and 2005 TN53. Values include 
the 1-cr uncertainty. (Epoch = JD2456200.5, 2012-Sep-30.0; J2000.0 ecliptic and equinox. Sources: JPL Small-Body Database and 
AstDyS-2.) 



2001 QR322 



2004 UP 



10 



2005 TN53 



semi-major axis, a (AU) 
eccentricity, e 
inclination, j (°) 
longitude of the ascending node, CI (°) 
argument of perihelion, a) (°) 

mean anomaly, M (°) 
perihelion distance, q (AU) 
absolute magnitude, H (mag) 



3().356±0.007 
0.02891±0.00010 
1.3210±0.0006 
151.548±0.015 
167.0±0.4 
60.3±0.4 
29.478±0.005 
7.7 



3().28±().()2 
0.03140±0.00010 
1.43107±0.00005 
34.798±0.003 
4±3 
345±2 
29.329±().()14 



30.262±0.()12 
0.069±0.002 
24.978±0.003 
9.2809±0.0002 
84.7±0.4 
296.7±0.2 
28.176±0.010 
9.1 



Table 3. Heliocentric Keplerian orbital elements of Neptune's L 4 Trojans (II): 2005 TO74, 2006 RJ103, and 2007 VL 305 . Values in- 
clude the 1-cr uncertainty. (Epoch = JD2456200.5, 2012-Sep-30.0; J2000.0 ecliptic and equinox. Sources: JPL Small-Body Database 
and AstDyS-2.) 





2005 TO74 


2006 RJ103 


2007 VL 305 


a (AU) 


30.262±0.013 


30.18 1 ±0.009 


30.195±0.014 


e - 


0.0536±0.0009 


0.0300±0.0010 


().0684±0.0003 


in 


5.248±0.002 


8.1615±0.0003 


28.1301±0.0010 


no 


169.370±0.006 


120.912±0.010 


188.5826±0.0010 


co(°) 


300.5±0.4 


17±1 


217.6±0.3 


MO 


278.7±0.4 


257±1 


359.2±0.2 


q (AU) 


28.639±0.011 


29.276±0.008 


28.130±0.012 


H (mag) 


8.5 


7.5 


8.0 



Table 4. Heliocentric Keplerian orbital elements of Neptune's L5 Trojans: 2004 KVig, 2008 LCis, and 201 1 HMio2- Values include 
the 1-cr uncertainty. (Epoch = JD2456200.5, 2012-Sep-30.0; J2000.0 ecliptic and equinox. Sources: JPL Small-Body Database and 
AstDyS-2.) 





2004 KVi 8 


2008 LC 18 


2011 HM102 


a (AU) 


30.108±0.011 


29.90±0.03 


30.05±0.05 


e - 


0.1847±0.0011 


0.086±0.003 


0.0785±0.0002 


if) 


13.6127±0.0015 


27.592±0.006 


29.4203±0.0004 


no 


235.6308±0.0005 


88.5219±0.0012 


100.98900±0.00006 


con 


294.2±0.3 


8±13 


152.05±0.03 


MO 


61.19+0.13 


173±15 


22.48±0.03 


q (AU) 


24.548±0.007 


27.327±0.007 


27.691±0.004 


H (mag) 


8.9 


8.4 


8.1 
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Fig. 4. Evolution of various parameters during the time interval (-500, 500) kyr for three Neptune L4 Trojans: 2001 QR322, 2004 
UP10, and 2005 TN53. The distance to the objects from Neptune is given in panel A; the value of the Hill sphere radius of Neptune, 
0.769 AU, is displayed. The resonant angle, A r , is displayed in panel B for the nominal orbit (source: JPL Small-Body Database 
and AstDyS-2). The orbital elements are depicted in panel C, a with the current value of Neptune's semi-major axis, panel D, e, 
and panel E, /. 2005 TN53 appears to be the most stable of Neptune's Trojans and 2001 QR322 the one with the largest libration 
amplitude. 
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Fig.5. Same as Fig.|4]but for Neptune L4 Trojans: 2005 TO74, 2006 RJ103, and 2007 VL305. 2007 VL305 exhibits the second largest 
libration amplitude among L4 Trojans. 
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Fig. 6. Same as Figs.g]and|5]but for Neptune L 5 Trojans: 2004 KVi 8 , 2008 LQs, and 201 1 HM 102 . The dynamical evolution of 2004 
KVi8 is clearly different and it is the most unstable of the "classical" Neptune Trojans. The most recent discovery, 201 1 HM102, is 
the most stable L5 Trojan and it may be as stable as 2005 TN53, likely the most stable of Neptune's Trojans. 
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